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LONGITUDINAL MASS-SPRING MODELING 
OF LAUNCH VEHICLES 

SUMMARY 


To determine longitudinal axisymmetric vibration modes of launch 
vehicles, simplified models must be established. As a first step, a launch 
vehicle can be idealized by a chain of masses, helical springs, and liquid-filled 
containers. However, the realistic mass-spring representation of the liquid 
propellants supported by the elastic container wall is the problem. In this 
report, application of Galerkin's method to the equations of motion of the above 
simplified system is proposed, using the container-liquid modes with special 
normalization factors as coordinate functions. This method, identical with the 
Ritz method in this case, is favorable in several aspects. First, the resulting 
model is based on the energy law. The use of the container-liquid modes, how- 
ever, results in accurate representation of significant structural and inertial 
characteristics of the container-liquid system. The special normalization of 
the coordinate functions converts the container-liquid systems directly in 
lumped parameter systems, having a number of degrees of freedom equal to 
the number of coordinate functions used. In this report, the containers are 
assumed cylindrical membranes with flat, rigid bottoms. For simplicity, 
typical tandem parts of vehicles, consisting of only one or two solid masses 
and springs and one liquid-filled container, are examined. If one-term 
approaches are considered, the resulting model can be interpreted as an actual 
spring-mass model. Because of the rapid decrease of the equivalent modal 
masses (with increasing mode number) , restriction to one-term approaches 
does not mean considerable loss of accuracy. Comparison of analytical results 
with experiments provides satisfactory agreement. In this way the use of large 
matrices in the final eigenvalue problem can be avoided. 


INTRODUCTION 


The "POGO" oscillations of launch vehicles are of primary interest. 
POGO oscillations result from a closed-loop interaction between one of the 
axisymmetric longitudinal vibration modes and the propellant system. Hence, 



the POGO problem has led to increased interest in the accurate determination 
of longitudinal axisymmetric modes and frequencies of vehicles. However, 
knowledge of the elastic and inertial characteristics of the structural components 
is required before the modes and frequencies can be calculated. Because of the 
complexity of launch vehicles , adequate simplified mass and stiffness models 
must be established for calculation of the desired modes. Usually only a few of 
the fundamental modes are of interest; thus, the structure can be idealized by 
taking only significant dynamic properties into account. In general, the 
representation of the solid masses is not as critical as that of the liquid pro- 
pellants. From experience it is known that the solid masses can be lumped at 
certain points of the vehicle axis and that the elasticity of the structure can be 
idealized by helical springs between the mass points. In this way, the solid 
structure is simplified by a chain of springs and masses. The liquid masses 
inside the elastic container must be treated differently. The liquid propellant 
masses constitute a high percentage of the overall vehicle masses throughout 
much of the power flight and, coupled with the solid vehicle masses and springs, 
may generate the fundamental modes of the complete vehicle. Thus , the correct 
comprehension of the liquid behavior inside the elastic container is of extreme 
importance. There is, however, no direct way to replace the liquid and the 
elastic container wall by lumped solid masses and springs. The problem of 
interaction between liquid and elastic container wall, a so-called hydroelastic 
problem, is governed by the laws off luid dynamics and elasticity. Its ana lytical 
formulation represents a boundary value problem. To ensure proper_represen - 
tation and clear understanding of the liquid-container interactions, thisjproblem 
must be solved. 


Nevertheless, one feels intuitively that the liquid masses supported by 
the elastic walls act like spring-mass systems coupled to the solid vehicle 
structure. Thus, it is understandable that engineers like to simulate the dynamic 
behavior of the liquid inside the elastic tanks by single spring-mass systems. 
Models of this kind have existed for a long time. Some of them are described 
in publications [1-5]. These models are based on the following assumptions: 

1. The equivalent vibrating solid mass equals or can be larger [3] 
than the total liquid mass. 

2. The liquid velocity field inside the tank, generated by the accelerated 
container motion, is neglected. 

3. Bending stresses and inertia of the wall are neglected. 
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In Appendix A a short description of these models is given, and their 
usefulness is shown by comparison with test results. 


Another kind of approach is presented in Archer [6] and Rubin [7]. 

Based on orthotropic shell theory, a finite element technique is utilized to 
construct the total launch vehicle stiffness and mass matrices by dividing the 
structure into ( axisymmetric) shell components, fluid components, and spring- 
mass components. The stiffness and mass matrices for the complete vehicle 
are obtained by superposition of the stiffness and mass matrices of the individual 
components, which are computed using a Rayleigh-Ritz approach. 


This paper is concerned with the modal analysis of vehicles , based on 
spring-mass modeling. However, the method proposed is quite different from 
that noted at the beginning. Analysis and model are based on the equations of 
motion of the reduced vehicle, consisting of a chain of masses, springs, and 
liquid-filled containers. Application of Galerkin's method reduces these equa- 
tions to a system of finite ordinary differential equations that can be interpreted 
as the equations of motion of a lumped parameter system, provided properly 
selected coordinate functions are used. This is explained briefly in the following 
paragraph. 


The equations of motion contain a finite number of ordinary differential 
equations (describing the motion of the solid mass-spring system) and a number 
of partial differential equations governing the liquid-filled container. The latter 
represent distributed systems having infinite degrees of freedom. To replace 
these equations by a number of ordinary differential equations, Galerkin's 
technique using the vibration modes of the partially liquid-filled container as 
coordinate functions is applied. In doing so, one obtains a finite number of 
coupled ordinary differential equations for each container and, collectively, a 
finite number of equations of motion of the entire vehicle. These equations can 
be written in matrix form displaying stiffness and mass matrix. In this way, 
the problem of free vibrations is reduced to a matrix eigenvalue problem. 


In the following sections it is shown that the container-liquid equations 
can be considered the equations of motion of lumped mass systems , provided 
the coordinate functions are normalized by specific factors. If this is the case. 


Galerkin's approach converts the container-liquid systems directly into equiva- 
lent lumped mass systems, whereby the number of masses (degrees of freedom) 
equals the chosen number of coordinate functions. In this way, relations evolve 
between the masses of the solid structure and those masses representing the 
liquid. 


Galerkin's method using the container-liquid modes as coordinate 
functions is also favorable from other aspects. Galerkin's approach, identical 
with the Ritz method in this case, ensures proper adherence to the law of energy, 
while the use of the modes results in accurate representation of significant 
structural and inertial properties of the container-liquid systems. The entire 
vehicle is then represented by a lumped parameter model, established without 
any assumptions violating the laws of fluid dynamics and elasticity. 


Thus, studying this model, one obtains understanding of the vibration 
problem at hand. First, one recognizes that the resulting model, in general, 
does not represent an actual spring-mass model. If interpreted in this way it 
would display "springs" having negative constants. However, the point is, 
if one-term approaches are considered, the container-liquid models represent 
actual spring-mass systems which agree, as far as the arrangement of masses 
and springs are concerned, with the model mentioned at the beginning, although 
the amount of mass and the spring constants are different. In this way, the 
entire vehicle is represented by a lumped spring-mass model. As shown 
below, the restriction to one-term approaches does not mean considerable loss 
in accuracy. One fact should be emphasized. To apply the approach suggested, 
the container-liquid modes must be known. This seems to be a disadvantage 
of the approach. However, in general, one cannot expect to be able to h andle 
the vibration problem of a complex structure without being able to handle the 
vibration problems of its single parts . 


To realize the model outlined above, specialization of the coordinate 
functions by normalizing factors is necessary. This procedure is based on the 
well-known principle of "Mechanical Analogy. " From the vibration modes, 
simple mechanical one-degree-of-freedom systems can be derived that respond 
in the same fashion as the modes , having the same energy and producing the 
same support forces. In Graham and Rodriguez [8] and Schmitt [9, 10] , this 
principle is applied to the lateral sloshing motion of liquid in a rigid container. 
Because of the special nature of this case, where the potential energy stems 
from the liquid surface elevation, a pendulum analogy results. Pozhalostin [11] 
contains a discussion of this principle, based on the longitudinal vibrations of 
an elastic container partially filled with liquid. In all these cases, forced 
vibrations were considered to derive the analog mechanical systems. 
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Thus, the single container-liquid modes must be studied. From the 
above remarks it can be concluded that any axisymmetric mode of vibration is 
equivalent to a single spring-mass system resting on the container bottom. The 
amount of mass and spring constant are determined by the following require- 
ments: potential and kinetic energy of the spring-mass system equal, at any 
time, those of the container-liquid system, vibrating in that particular mode. 

In addition, the resulting liquid pressure force on the container bottom and the 
force exerted there by the mechanical spring are equal. Consequently, the 
particular mode used must be normalized to correspond with the container 
geometry and liquid height. In the following section, the equivalent single 
spring-mass systems of the modes of a partially liquid-filled cylindrical 
membrane with flat, rigid bottom are derived. It is easy to confirm that the 
sum of all "modal masses" determined as outlined above is smaller than the 
total liquid mass. Furthermore, one realizes the rapid decrease of the modal 
masses with increasing mode number. 

Using these normalized modes as coordinate functions, the Galerkin 
approach can be performed. Because of the rapid decrease of the modal 
masses, only a few degrees of freedom must be added to obtain a realistic 
model of the container-liquid system. In this way, the use of large matrices 
in the final eigenvalue problem can be avoided. In the third section, the equa- 
tions of motion of typical tandem parts of vehicles consisting of solid masses, 
springs, and a partially liquid-filled cylindrical membrane are derived. In the 
following sections, Galerkin’s method is applied to two of these cases: partially 
liquid-filled cylindrical membrane with a top mass; partially liquid-filled 
membrane resting on a single spring-mass system. Two special cases were 
numerically solved and compared with test-results performed by Southwest 
Research Institute, San Antonio, Texas, under Contracts NAS8-11045 and 
NAS8- 20329. As seen in the last two sections, the agreement is satisfactory. 

The solutions of these special cases render possible the spring-mass modeling 
of any chain of solid springs and masses and liquid-filled cylindrical membranes. 


EQUIVALENT SINGLE SPRING-MASS MODELS 
OF THE VIBRATION MODES 


The axisymmetric free vibrations of partially liquid-filled cylindrical 
membranes with flat, rigid bottoms are considered. The membrane is referred 
to a cylindrical system of coordinates r, 6 , z (Fig. la). 
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a) Coordinate System; Displacements b) Equivalent Single Spring- 

th 

Mass System of the j Mode, 

FIGURE i. PARTIALLY LIQUID-FILLED CYLINDRICAL MEMBRANE 
The following assumptions are made: 

1. The wall inertia is neglected. 

2. The membrane is free at the top. 

3. The liquid is considered inviscid and the flow irrotational. 

4. The pressure created by the surface wave is neglected. 

From assumptions i and 2, the equation of motion of the membrane 
[13-15] follows as 

w(z, t) = p(R, z, t) . (1) 


From assumption 3, it follows that the liquid velocity inside the tank can 
be derived from a potential <p( r, z, t) , which is the solution of Laplace equation. 


V 2 c p = 


9 2 <fr 

3r 2 


1 dcp 9 2 0 
r 9r + 9z 2 


0 < r < R 
0 < z < h 


( 2 ) 
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The liquid pressure caused by the velocity potential is then given by 


p(r, z, t) 


. 80( r, z, t) 


( 3 ) 


The boundary conditions are 


z = 0 


r = R 


80 

8z 


= 0 


80 _ 8w 
8r = 8t 


z=h , 0=0 


(4) 

(5) 

( 6 ) 


Condition ( 6) follows from assumption 4 and equation ( 3) . 

As well-known [ 14, 15] , equation ( 2) and the boundary conditions ( 4) 
through ( 6) represent a boundary value problem with the following solution: 


0X r, z, t) = a^(t) $At,z) 


V r,Z)= °j R Io(X j h> C ° s(X j 
(2j-l)ir 

2 »j — 4,2, 3, . . . 


(7) 


( Proof by substitution. ) The quantity a^( t) has the dimension of a velocity 
and is determined by initial conditions. The dimensionless constant, c., will 

J 

be determined in the following. According to equations ( 3) and ( 7) the liquid 
pressure created by 0^( r,'z, t) follows as 


80. (r, z, t) 


Pj (r, Z ,t) «-p— 


= -pa ; (t) 0.(r, z) 
1 1 


( 8 ) 


It can be concluded easily that the solution (7) fulfills boundary conditions (4) 
and (6). 
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From boundary condition (5) and equation (7) the 
displacement follows as 


w.(z, t) = a.(t) w.(z) 
3 J 3 


w.(z) 

3 


9<jMR, z) 

J — — 

9r 


c. A. — I^A.-jr ) cos (A. 

J J h 1 j h j 

j = 1> 2, 3, ... 



radial membrane 


(9) 


Substitution of equations ( 7) , ( 8), and ( 9) into equation ( i) yields 


( 10 ) 


(ID 


til 

Equations (7) , ( 9) , and ( 11) represent the j solution of the eigenvalue 
problem at hand: the free vibration of the partially liquid-filled cylindrical 
membrane having a flat, rigid bottom. Since the problem is homogeneous, the 
solution is determined only up to the constant factor c. as shown by equations 

th 1 

(7) and ( 9) . As noted in the introduction, the j vibration mode is equivalent 
to a single spring-mass system resting on the container bottom and having, at 
any time, the same kinetic energy, potential energy, and reaction force on the 
container bottom. From the latter three conditions, the amount of mass and 
the spring constant of the mass-spring system and the modal factor c. can be 
determined. ^ 

Next, the law of conservation of energy of the liquid-filled membrane, 

vibrating in the mode, is examined. Multiplication of equation (10) by 
9w.(z, t)/9t (given by equation ( 9) ) integration over the wetted membrane wall 

and over a certain time interval gives 


E6 - 


P0.(R, z)a. + — ry w.(z)a. = 0 
J J R 2 J 


or 


_j_ 2 E6 

-*- = i 0. = '5 

a. j R S P 


1 


A. Ii(A. p) 

3 h J h _ 


R 

j h 


I„u. 


2R?r 

2 



£ (R, z) 


w.(z) a. 2 + 
3 3 


2Rt r E6 
2 R 2 


h 

f [w.( z)] 2 dz a. 2 
0 J J 


= const. 


( 12 ) 
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Equation ( 12) expresses the law of conservation of energy of the system 

til 

membrane-liquid, vibrating in the j mode. It can also be assumed to be the 
energy law of a single spring-mass system having the mass 


— R R R 

m. = 2Rirp J 0.(R, z)w.(z) dz = c } A. I<j(A. — > Ijl X. — )m T (13) 

J n j 3 33 3 3 


and the spring constant 


k. 

3 


■ 2R,r |f /Iw.(z)] 2 dz = “ [c. A.IjU. f)} 2 


2Rtt E6 
h 


(14) 


The equation of motion of the single spring mass system (Fig. lb) follows then 
from equations ( 12) 'through ( 14) as 


m. a. + k. a. = 0 ( 15) 

3 3 3 3 

Because it is assumed that the spring of this system rests on the container 
bottom and that the spring force equals the force exerted there by the liquid, it 
follows 


R 

-k. a. = 2 5r f p.(r,o, t)rdr 
3 3 -jj 3 

or, under consideration of equations (7) , ( 8) , and (15), ([16], p. 484) 


R - 2 R 

m^ = - 2irp f <p.( r, o) rdr = -c. - I^A. — ) 


m. 


( 16 ) 


Thus, from equations (13) and (16) one obtains 


V 


2 _ 

A. 




(17) 
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Now equations (7) , ( 9) , ( 13) , 


( 14) , and ( 17) result in 


r 2h I ° (A jh ) ,, z x 

^ oos<x ^ ) 


jh 


- 2 

W j "A. T/ , R, 

j i„u jir ) 


cos( A. £) 
J h 


4 I ' a i r ) 
m. = — 


j X? , K r(1 R, 


m r 


k.= 2 

3 






2Rtt E6 


j - If 2, 3, . . 


( 18 ) 


(19) 


( 20 ) 


( 21 ) 


Equations (15) , (20) and (21) define completely the equivalent single 

th 

spring-mass system of the j vibration mode. From equations ( 20) and (21) 
also follows the circular frequency as given by equation ( 11) . 

To determine the spring-mass model numerically, the function, 


Il(x) 

Io(x) 


, 0 ^ x ^ °° 


( 22 ) 


must be evaluated. As shown in Appendix A, continued fractions and asymptotic 
expansions prove to be useful for that purpose. 

Figure 2 is the modal mass ratio ( 20) , 


m j 

— , j = 1, 2, 3 
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versus h/R. It manifests the importance of the first mass mj and the rapid 
decrease of the higher modal masses. 


From equation (A-i) , it follows 


lim 

h-*-°° 




R 


) 


A j h Io(A j h ) 


1 

2 


(23) 


whereby the limes i/2 also represents the maximum value. Thus, from equation 
( 20) it can be concluded 


hence 


m. 


lim sup 


m r 


lim sup 


m„ 


lim sup 


m? 

m. 


lim sup 

m T 


8 

( 2j - 1) 2 7T 2 


= 0. 8106 
= 0.0901 
= 0.0324 


» 


The above values, also in Figure 2, demonstrate again the degree of 
significance of the single modal masses. 

A point of primary interest may now be touched upon. The amount of the 
th 

j modal mass is given by equation ( 20) . From this equation and the third 
equation ( 7) , the sum of all modal masses follows as 


OO 


2 

j=i 


m. = 
J 


16 

— r-m„ 
ifi T 



**} £> 


( 2 j - 1) 2 R R 

jh IoU jh } 


Now, the question arises whether 


(24) 


oo 
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From equations (23) and (24) it follows ([16], p. 808) 


lim 
h— 00 




CO 

— _JL V 1 
m j ~ it 2 ( 2j - l) 2 


1 


Because 1/2 represents the maximum value of the function in equation ( 23) , 
which never can be reached, one concludes 


E m , <m T 

j-i 1 

The following integrals should be mentioned for later use: 
i h 

— f w.( z) dz = (-1)- 1 k . 

R n 3 J 


R 


2 ^Jh 1 


j = X. . R _ - . R . 
j K T In( K T> 
J j h u j h 


2Rirp f (h - z) w.( z) dz = -m. 
n J J 


0 


where w/ z) and rq. are given by equations (19) and (20) . 


(25) 


(26) 


EQUATIONS OF MOTION OF TYPICAL TANDEM PARTS 

OF VEHICLES 


Figure 3 shows a tandem part of a vehicle consisting of two solid masses, 
springs , and a partially liquid-filled cylindrical membrane. 

To obtain the equations of motion, the system shall be divided into sub- 
systems for which the equations of dynamic equilibrium can easily be set up. 
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a) TANDEM PART OF A VEHICLE b) EQUIVALENT LUMPED PARAM- 
CONSISTING OF TWO SOLID ETER SYSTEM. 

MASSES, TWO SPRINGS, AND A 
PARTIALLY LIQUID- FILLED 
CONTAINER. 

FIGURE 3. 

In doing so, the membrane is crossed perpendicular to the axis somewhere 
between the masses M 2 . The equilibrium conditions of the lower and upper 
part are given by 


R 


Mjxj + Kjxj - 2 RttN + 2t r f p(r, o, t) rdr = 0 

(27) 

2 0 


M 2 x 2 + K 2 x 2 + 2R7T N = 0 , 

^ z 

(28) 


where 

Xt = Xj( t) , x 2 = x 2 (t) 
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If the wall inertia is neglected ( assumption 1) , the equilibrium condition of a 
membrane ring element one inch high is given by 



0 £ z < h 
h < z < 1 


(29) 


The membrane forces N , N„ can 
. ,, z 6 

as follows: 


N = D 
z 


V D 


/9u A w \ 
(az v R ) 

/ w 9u \ 
\R + v dz ) 


u =u( z, t) ; w = w( z, t) ; D 


E6 

i-v 2 


be replaced by the displacements [13] 



(30) 


Substitution of equation ( 30) into equations ( 27) through ( 29) results in the 
three equations of motion containing the four unknown quantities u, w, x lt x 2 . 


To reduce the number of unknowns, u shall be eliminated. Starting from 
equation ( 30) , the displacements can be expressed as follows: 


All 

E<5 ~ - = N - i/N 

dz z d 


E6 77 = N_ - i^N 
R 6 z 


(31) 


Taking into account the independence of N on z in <^0, z^> , as indicated by 
equations ( 27) and ( 28) , and the fact that N = 0 in <(h, f> , one obtains from 

U 

equations (31) , by integration over the intervals <b, iy and <o, h> , 
respectively, 

h 

E6(x 2 - Xl ) = lN z - v J N 0 dz 



7 “ dz = - vh N + 
R z 




dz 


» 
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and finally 


N 

z 


kp 

2Rtt 


( x 2 - Xj + v 



k 

o 


2Rtt E6 
Ul-V 2 j) 


( 32 ) 


From equations ( 27) through ( 29) , ( 31) ( second equation) , and ( 32) , the 
equations of motion ( now independent of u) follows as 


h R 

(K! + k 0 ) xj-^kp J ^ dz - k 0 x 2 =-M 1 x 1 - 2ir f p(r, o, t)rdr 
0 K 0 

h 

-j'k 0 x 1 + 2R7rE6 ^- + v 2 k 0 f “ dz + vk 0 x 2 = 2R 2 7rp(R, z, t) 

R o R 

h 

- koXi + t'kp f ^-dz + (K 2 + k 0 )x 2 = -M 2 x 2 


(33) 


where 


X! = xj(t) ; x 2 = x 2 (t) ; w = w(z,t) ; 0 ^ z < h 

Under assumptions 3 and 4 of the preceding section, the liquid flow is 
governed by a velocity potential <p(r, z, t) . In accordance with equation (3) , the 
pressure can be expressed by 


p(r, z, t) 


= -P 


90(r, z, t) 

at 


- (h - z) Xj 


(34) 


To solve equation (33) by Galerkin's approach, using the container- 
liquid modes of the preceding section, one has to assume, in accordance with 
equations ( 3) , (7) through ( 9) , ( 9 ) 
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( 35 ) 


n 

<M r, z, t) = Yj fA.(t) - k i (t)]<p.{r, z) 

j= i 3 3 


w(z, t) = Y lA.(t) -x^tjjw.tz) 

3 = 1 3 3 

where <£.( r, z, t) and w.(z) are given by equations ( 18) and ( 19) . Then the 
well-known procedure can be performed. 

Before Galerkin's method is applied to the system of equations ( 33) , two 
simpler systems will be examined. 



1) Figure 4 is a cylindrical membrane partially filled with liquid and 
fixed on the lower end having a top mass. 



a) PARTIALLY LIQUID-FILLED b) EQUIVALENT LUMPED PARAM- 
C YLINDRICA L MEMBRANE ETER SYSTEM 

HAVING A TOP MASS 


FIGURE 4. 
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The equations of motion follow from the second and third equations ( 33) 
by setting 


Xj = 0 

K 2 = 0 
One obtains 


2RirE6 ^ + I'kolxj+y J dz) = 2R 2 7rp(R, z, t) 


R 


ko(x 2 + v j — dz) = - M 2 x 2 


(36) 


(37) 


From equations ( 34) and ( 36) , it follows 




(38) 


2) Figure 5 is a partially liquid-filled cylindrical membrane resting 
on a single spring mass system. 

The equations of motion follow from the system of equation (33) by 

setting 

K 2 = 0 
M 2 = 0 
One obtains 


(39) 


R 

Kjxj = -MjX! -2ir Jp(r, o, t)rdr 
0 


W 

E<5~ = Rp( R, z, t) 


( 40 ) 


The liquid pressure is given by equation ( 34). 
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a) PARTIALLY LIQUID-FILLED b) EQUIVALENT LUMPED SPRING- 

CYLINDRICAL MEMBRANE MASS SYSTEM 

RESTING ON A SINGLE SPRING- 
MASS SYSTEM 


FIGURE 5. 


Equation ( 39) expresses the fact that the upper end of the membrane is 
free. Thus, no axial membrane forces are acting during vibration of the 
system. Consequently, equations (40) do not contain any terms stemming 
from Poisson effects. 

In the following, Galerkin's method using the coordinate functions 
of equations ( i 8) and (19) will be applied in equations ( 33) . From the result- 
ing linear system of ordinary differential equations, which represent the 
dynamic equilibrium conditions of the mass points involved, the lumped mass- 
spring model sketched by Figure 3b will be derived. 

Equations ( 34) and ( 35) must be substituted into equations ( 33) , and 
the specified integrations must be carried out. Then the second equation ( 33) 
must be multiplied by w.( z) , (i = i, 2 . . n) and integrated over <o, h> ; 

whereby equations ( 13) , ( 14) , ( 16) , ( 25) , and ( 26) must be taken into account. 
Diagonalizing of the mass matrix requires simple intermediate operations, 
which are not shown. One obtains finally 


I 
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and 


ko 

are given by equations ( 20) , ( 21) , ( 25) , and ( 32) . 

Furthermore , 

k* = 1 + v Yj (-1)^*. • 

5 = 1 3 

Because the mass matrix of equation (41) is diagonal, any of the 
equations represents the dynamic equilibrium condition of the particular point 
mass appearing in that equation. The stiffness matrix contains no zero element; 
thus, one concludes that all masses are connected with one another by springs. 
Figure 3b is a spring-mass system having such an arrangement of masses and 
springs. Massless rigid bars drawn in Figure 3b serve for connection of all 
masses with one another by the vertical springs. The stiffness matrix of this 
system can easily be determined as 
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Comparison of the elements of this matrix with those of the stiffness 
matrix of equation (41) results in a system of linear equations having the 
unique solution 


= (-1)^ 1 vko fc 


K.. = K. . = ( -1) 
K (0) = k 0 K* 


i+i-1 2i 

V knK.K. 

1 J 


= k. + (-l)-* vVl 0 k*k. 


i*J 


ijj 1, 2, ... n 


(43) 


From these equations it can be seen that the "spring" constants and 
K.. are negative if j and i + j, respectively, are even numbers. The "springs" 

referring to these constants are dotted in Figure 3b. Consequently, the lumped 
system of Figure 3b represents an actual spring -mass system only if j = 1 
( one-term Galerkin approach) . Because of the predominance of the first modal 
mass, restriction to the one-term approach will be sufficient in most cases. 


As already mentioned, the arrangement of masses and springs resulting 
from the one-term Galerkin approach is the same as that of Wood's model, 
(Appendix B) . Of course, the amount of the equivalent liquid mass and the 
spring constant differ considerably. 


PART I ALLY LIQUID-FILLED CYLINDRICAL MEMB RANE 
HAVING AMASS ON TOP 


The free vibrations of this system, shown in Figure 4a, are governed 
by equations (37) and (38). As in the previous case, application of Galerkin's 
method results in a system of ordinary differential equations of motion having 
stiffness and mass matrix as 
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k 0 -t'koKj (-D^ko^ 

-t'koKt ki+i^koKi (-l) i + r V 2 k 0 K 1 K Q 

Vk 0 K 2 -V 2 VL 0 K 2 K i (-l) 2+ Vk 0 K 2 K n 

(-i) n ^k 0 K (-i) n+ 1 v 2 k 0 k Ki k +^ 2 k 0 (c 2 

n. n n n 

M 2 0 0 o 

0 mi 0 0 

0 0 m 2 0 


0 0 0 m 

n 


This result can also be obtained directly from equations ( 36) and ( 4i) . 
The equivalent lumped parameter system is represented by Figure 4b. 

Because of the sharp decrease of the higher modal masses, as shown 
in Figure 2, their coupling effect is small. Thus, if higher modes of the 
coupled system are not asked for , consideration of only the first modal mass 
will, in many cases, result in sufficient accuracy. Then equations (41) reduce 
to the simple system 


k 0 -vk 0 K i 


*1 

-vkflKi kj + v\ 0 k\ 


_ a L 


1 ~~ 

O 

CM 

l5_ 


Xl 

L° m i_ 


_a L 


(44) 


having a quadratic frequency equation. From equations ( 42) and ( 43) , the 
spring constants of this model follow as 
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9 


(45) 


K (0) = k 0 (i - VK t ) 

K (l) = ^k 0 Ki 

K (l) = kt - vk 0 Ki + v 2 ko K i 

where k ls K t , and k 0 are given by equations ( 21) , ( 25) , and ( 32) . Equation ( 20) 
gives m*. 

To prove the accuracy of the one-term approach, analytical and experi- 
mental results using a stainless steel container with three different top masses 
M 2 (weight W) were compared. The results are shown in Figures 6 through 8. 
The agreement is satisfactory. 


PARTIALLY LIQUID-FILLED MEMBRANE RESTING 
ON A SINGLE SPRING-MASS SYSTEM 


Figure 5a is the vibrating system which is governed by equations (34) 
and ( 40) . The linear system of ordinary differential equations of motion 
follows by application of Galerkin's technique, or it can be concluded directly 
from equations ( 39) and (41) . In the latter case, one has to consider that all 
terms containing v vanish. 


Stiffness and mass matrix follow as 


ki 

0 .... 

0 

~kl 

0 

k 2 .... 

0 


0 

0 .... 

k 

n 

-k 

n 

a 

■ki 

-k .... 
n 

-k 

n 

K i + 2 

j-i 


mj 0 


0 


m 2 0 


0 


0 

0 


0 m 


n 


0 0 M- 


n 

+ m - £ m 

J-l U 
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Natural Frequency, cps 


-| — Ft H“ Model of the 

1-term approach 

K = 5.0 in. (12.70 cm) 

l = 15.0 in. (38. 10 cm) 

6 = 0.005 in (0.0127 cm) 

E = 30 x 10 6 lb /in . 2 ( 2. 100 x 10 11 N/m 2 ) 
p = 9.359 x 10" 5 Ib-sec 2 /in. 4 
(10 3 kg/m 3 ) 


\ T 


Second Mode 



1 — i 


First Mode 


— Experiment 

— 1-term Galerkin Approach 

— Wood's Model W 2 


JL 


Liquid Depth Ratio, h/1 


FIGURE 6. NATURAL FREQUENCIES OF PARTIALLY FILLED TANK WITH 

TOP MASS W = 11. 48 lb ( 51. 066 N) 





FIGURE 7. NATURAL FREQUENCIES OF PARTIALLY FILLED TANK WITH 

TOP MASS W = 22. 32 lb ( 99. 284 N) 
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Natural Frequency, cps 


FIGURE 



8. NATURAL FREQUENCIES OF PARTIALLY FILLED TANK WITH 
TOP MASS W = 34. 53 lb ( 153. 597 N) 



where m. and k.(j = 1, 2, . . . n) are given by equations ( 20) and ( 21) . 

3 J 

Because all Poisson terms are zero, it follows from equations (42) and 

(43) 


= Ky = 0 i, j = 1, 2 ... n 

= k. j = 1, 2, 3, . . . n 

3 


In this case, no "springs" with negative constants exist. Consequently, 
the equivalent system is an actual spring-mass system for all values of n. It 
is represented by Figure 5b. 


Because the configuration of the nonzero elements of the stiffness matrix 
is relatively simple, the frequency equation can be set up immediately. For a 
one-term and a two-term approach these equations are 


X 2 - 


X 3 - 


2 , 2 
CJq + Cdj 


2 2 

X + <jl>q co 1 = 0 


4 + u? ( 1 + a») +6 j( 1 + 5 ,) 

2 2 , 2 2 /. . mim2 \ 2 2 

O>0 wi + Wi W 2 + -jj" + +U 0 U2 


, 2 2 2 n 
X - coo wj co 2 = 0 


where 


CO 


2 

0 


Ki 

M 


M = Mj + m^ 



m. 

J 


and 


n = 1, 2 


(46) 


X - ( 2rrf) 2 


(47) 
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Table I is the test and analysis results from a one -term and a two-term 
approach at different filling heights. The first frequency resulting from the one- 
term approach is already as good as that of the two-term approach. Test values 
of the third frequency could be obtained for the highest levels only. This fact is 
not surprising, however, in view of the smallness of the second modal mass. 

The container used is a ring-stiffened Mylar cylinder with rigid bottom. The 
effective modulus of elasticity was determined by use of a correction factor as 
outlined in [ 12] . 


CONCLUSIONS 


Longitudinal spring-mass models of launch vehicles can be obtained 
by application of Galerkin's technique. As a first step, the vehicle can be 
idealized by lumping its solid masses and (longitudinal) springs to a chain of 
lumped masses, springs, and liquid-filled containers. The equations of motion 
of this vehicle model consist, then, of a number of ordinary differential equations 
describing the motion of the solid point masses, coupled with partial differential 
equations governing the liquid motion inside the container. Application of 
Galerkin's technique using the container-liquid modes with special normalization 
factors reduces this system of equations to the equations of motion of a finite 
degree of a freedom -lumped parameter system. 

In this report, the containers are assumed cylindrical membranes 
with flat, rigid bottoms. For simplicity, typical tandem parts of vehicles 
consisting of only one or two solid masses, springs, and one liquid-filled 
container are examined. From the analysis the following can be concluded: 

1. The lumped parameter systems resulting from one-term approaches 
can be interpreted as actual spring-mass systems. Thereby the container- 
liquid system is replaced by one modal mass and three helical springs ( Figs. 

6 through 8) . Spring constants and modal mass are given by equations ( 20) , 

( 21) , ( 25) , ( 32) , and ( 45) . Because of the rapid decrease of the equivalent 
modal masses, no considerable loss of accuracy results from one-term 
approaches. This behavior is proven by comparison of analytical results with 
experiments (Figs. 6 through 8, Table 1). 

2. If approaches with more than one term are performed, the resulting 
lumped parameter system can no longer be considered an actual spring-mass 
model. In this case, the spring-mass arrangement that replaces the container- 
liquid subsystem exhibits "springs" with negative constants descending from 
Poisson effects. These springs are sketched with dotted lines in Figures 3 and 
4. Spring constants and modal masses are given by equations ( 20) , ( 21) , ( 25) , 

( 32) , ( 42) , and ( 43) . 
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TABLE I. NATURAL FREQUENCIES OF A PARTIALLY LIQUID- FILLED 
MYLAR CONTAINER RESTING ON A SINGLE SPRING-MASS SYSTEM. 
ONE-TERM APPROACH; TWO-TERM APPROACH; TEST VALUES. 


= 11. 2 in. (28.45 cm) 

= 4. 5 in. ( 11. 43 cm) 

. lb sec 2 
= 9. 359 x 10 s — 7^-4 


K = 11. 62 — (2. 035 X 10 e — 1 
in. m 


W =4.67 lb (20.773 N) 

E 6 = 5050 7 ^ (884390 — ) 
eff in. m 


R -i 



Mj + m - m t 


(— ! m 2 

T k i f k r 

M 1 +m T -m 1 -rn 2 


one-term appioach model two-term apptoach model 


h 


ft 



*2 



h 



1-Term 

2-Term 

Test 

i-Term 

2 -Term 

Test 

1 

1-Term | 

2-Term 

Test 

ii. i 

42. 5 

42. 5 

43. 1 

131. 8 

115. 4 

116. 0 


173. 9 

167. ( 

10. 1 

45. 6 

45. 6 

46.2 

136. 6 

122. 3 

122. 0 


180. 7 

175. ( 

9. 1 

49.2 

49. 2 

49. 6 

142. 1 

129. 8 

129. 0 


188.7 

- 

8. 1 

53. 4 

53. 4 

53. 9 

148. 4 

138. 3 

136. 0 


198. 2 

194. ( 

7. 1 

58. 4 

58. 3 

58.7 

154. 5 

147. 0 

140. 0 


209. 4 

_ 






Although this paper deals only with cylindrical membranes having a flat, 
rigid bottom , the method outlined can be extended to general axisymmetric 
membrane or shell containers, provided the container-liquid modes are known. 

In reality, many tanks consist of long cylindrical parts and spherical or elliptical 
ends. During much of the power flight time, even if the propellant masses are 
still large, the container-liquid system will behave like a liquid-filled cylinder 
with a flat, rigid bottom. Certainly an elliptical or spherical end represents a 
disturbance of the liquid flow, and its influence must be investigated. However, 
this is only one of many other disturbances present in the container-liquid sys- 
tem. Especially the behavior of tanks stiffened by rings and stringers, and also 
the influence of anti-slosh devices, must be studied. For cylindrical tanks with 
stiffening rings only, a correction factor of Young’s modulus can be derived so 
that the ring-stiffened wall is approximated by a wall of uniform thickness having 
a slightly greater modulus of elasticity [12], 
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APPENDIX A 


NUMERICAL EVALUATION OF Ijlxl/lglx) 


In the section on equivalent spring-mass model of vibration modes , 
normalized modes and equivalent masses and springs are represented in terms 
of I^xJ/IqCx) . For numerical evaluation, useful expressions of this function 
should be considered. 

For argument values 

x < 10 

a continued fraction representation proves to be convenient. In [16] the follow- 
ing expansion is given: 


= L_ 

Jo<«> 2 _JL 



iL 

8 - . 


where Jo( £) and J*( |) are the Bessel functions of orders 0 and 1 of the first 
kind. 


Setting £ 2 = -x 2 , 

the above continued fraction transforms into 

Il(x) = x_ 

4(x) 9 , *1 

z+ 4 

( [16] p. 375 9. 6. 3) 



( A-l) 
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For argument values 


x — 10 


asymptotic expansions of the modified Bessel functions can be used [16]. 


ll(x) 
I 0 ( x) 


4-1 (4-1) (4-9) (4-1) (4-9) (4-25) 

_ _8x + 21 ( 8x) 2 3'.{8x) 3 _ 

7 "l , 1. 9" l-~9. 25 

+ 8x + 2!( 8x) 2 3J( 8x) 3 


Figure A- 1 shows I*( x)/I<)( x) , and 


Il( x) 
xl 0 ( x) 


versus x. 



FIGURE A-i. Ii (x)/Io (x) AND Ij (x)/xIo (x) VERSUS x 



APPENDIX B 


WOOD'S MODEL 


Figure B-l is Wood's model of a cylindrical tank with top mass and 
rigid bottom. As already mentioned, the arrangement of springs and masses 
agrees with that of a one-term Galerkin approach, while the amount of equiva- 
lent mass and the spring constants are different. There is, however, no unique 
interpretation of Wood's model. In addition to the original model, two modifi- 
cations can be found in the literature on the subject. Spring constants and 
equivalent masses of these models, designated by W lt W 2 , and W 3 , are gathered 
in Table B-I. Wood's model is not restricted to cylindrical tanks with rigid 
bottoms. Versions for tanks with flexible bulkheads and tanks stiffened by rings 
and stringers also exist [1, 4]. 

Obviously, important structural and inertial properties of the container- 
liquid system cannot be represented by simple models of this kind. There are, 
however, some other weaknesses of these models. Consider for instance the 
spfing constants of the models Wj and W 3 . If h goes to zero, all three spring 

constants K^i (j = 1,2,3) go to infinite, thus both frequencies of the models 
w 

go to infinite. This behavior is not in agreement with test results, as shown by 
Figures 6 through 8. Now consider the model W 2 . If M vanishes, the container- 
liquid frequency, determined by K w ^ and nij, goes to zero if the container 
length 1 goes to infinite. This is, however, inconsistent with theory 
and experience. 

In Figures 6 through 8, analytical results obtained by application of 
Wood's model W 2 are also presented and can be compared with experimental 
results. 
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a) PARTIALLY LIQUID FILLED b) WOOD'S MODEL 
CYLINDRICAL TANK WITH TOP 
MASS 

FIGURE B-l. 

TABLE B-I. SPRING CONSTANTS AND MASS OF WOOD'S MODEL 


Original Model 

Wi 


Modified Models 


2(1-!/) 

k 

3(1-!/) 

3-21/2 V 

4 - 3v 2 ^ 

4- 3i/ 2 

2p , 

3^ . 

3* y 

GO 

1 

to 

to 

4 - 3i/ 2 W 

4-3v2 11 

3 - 2i/ 

4- 3v 

4 - 3i/ , 

3-21/2 V 

4-3!/ 2 W 

4-31/2 k 

2Rtt E5 

2Rtt E6 

2R7T E6 

h 

1 

h 

R 2 7r hp 

> R 2 7T ph 

R 2 7t ph 

[1] [2] 

[3] 

[5] 
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